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' —,' ' Abstract 

^ ; 

Qh| The cross section for Higgs+jct production in bottom quark annihilation is calcu- 

O ■ lated through NLO QCD. The five-flavour scheme is used to derive this contribution 

D I to the Higgs+jet production cross section which becomes numerically important in 

y^/ the MSSM for large values of tan/3. We present numerical results for a proton collider 

with 14TeV center-of-mass energy. The NLO matrix elements for da/dpT are then 

combined with the total inclusive cross section in order to derive the integrated cross 

section with a maximum cut on px at next-to-next-to-leading order. 

^H ! 

in ■ 1 Introduction 

O 

o 



The Higgs mechanism [IH3] plays a central role in both the Standard Model (SM) [IHB] 
and its supersymmetric extensions [?]• The gauge bosons and quarks acquire masses 
through interactions with Higgs fields. Up to now, the search for the physical Higgs boson 
has been unsuccessful but has led to the exclusion of a certain Higgs mass range [8l[9]. 
In combination with the fits of electro-weak precision data to higher order perturbative 
5^ \ calculations this leads to a rather small range of allowed values for a SM Higgs boson 

mass [To]. 



Supersymmetric (SUSY) theories require an enlarged Higgs sector. The minimal SUSY 
extension of the SM leads to five physical Higgs bosons. Due to the larger number of free 
parameters in SUSY, the sensitivity of experimental data to Higgs bosons is weaker than 
intheSM \Tl[\T2\. 



The Large Hadron Collider (LHC) is expected to find a Higgs boson if it exists. To do this, 
various production and decay channels must be considered. The relative utility of each 
channel depends strongly on the Higgs mass and couplings. While in the SM gluon fusion 
is the Higgs production process with the largest cross section by far, in SUSY theories with 
large tan/3, Higgs production in association with bottom quarks is dominant (for reviews, 
see Refs. |131I14|: detailed studies of the relative importance of gluon fusion and bottom 
annihilation have been performed in Ref. |15fll7j ). This is because in this region of the 
SUSY parameter space the Hbb coupling is enhanced relative to the SM. 

Assuming that only four quark flavours and the gluon make up the proton (the so-called 
"four-flavour scheme" or 4FS), the dominant leading order Feynman diagram for this 



(a) 

Figure 1: Leading order diagrams for associated bhh production in the (a) four and (b) 
five flavour scheme. 




process is shown in Fig. 1(a) If the bottom quark was massless, integration over phase 
space would lead to divergences arising from the kinematical region where one or both 
bottom quarks are collinear to the incoming partons. The bottom quark mass regulates 
these divergences, but they still leave traces in terms of logarithms of the form ln(?n,^/m^). 
Such logarithms lead to large perturbative coefficients, so ideally one would like to resum 
them. This can be achieved by considering this process in the five-flavour scheme (5FS) [181 
[19], i.e. by introducing bottom quark PDFs (parton density functions) HI Now that the 
h quarks can appear in the initial state, the leading order process is changed to that 



in Fig. 1(b) We note that the scheme choice amounts merely to a re-ordering of the 
perturbative series. Of course, when truncated at a finite order, results obtained in either 
scheme will differ, with the difference being formally of higher order in as- 

However, it was found that the difference between the inclusive cross section in the 4FS 
and the 5FS differs by roughly a factor of five when evaluated at /ip = A*R = AfH) where 
/xf//^r is the factorization/renormalization scale. This remains true also at NLO QCD 
which was calculated for the 5FS in Ref. [2ll[22], and for the 4FS in Ref. [23l[24]. It was 
thus proposed in Refs. [22 1I25H27] that when using the five flavour scheme the appropriate 
scale choice is mff/4. 

A weakness of the 5FS is that it neglects the contribution from large-pr bottom quarks at 
leading order. This is taken into account only at NNLO and higher (note that the LO set of 
Feynman diagrams in the 4FS is part of the NNLO set in the 5FS). Indeed, the factorization 
scale dependence at NNLO is very flat [28] and seems to confirm the "natural" scale choice 
at lower orders of /ip = -^h/4. 

It has been pointed out a long time ago that it can be advantageous to consider the i?+jet 
process instead of the fully inclusive production when searching for the Higgs boson [29] . 
The gg —^ H+^et cross section is known at LO, including the full top and bottom quark 
mass dependence, both in the SM [29H3T] and in the MSSM [31]. NLO QCD corrections are 
known in the heavy-top limit |32H34] . It is expected that a very good approximation of 
the MSSM effects can be obtained by simply replacing the corresponding Wilson coefficient 
of the effective ggh coupling with its MSSM expression [35H37] . at least as long as tan/3 
is not too large. Otherwise, bottom loop effects which are not covered in the heavy-top 
limit will be important. Resummation effects for small and large px of the Higgs boson 
have been treated in Refs. [SHHiTJ. 



As mentioned already above, for large tan/3, it is essential to also take bottom quark 
annihilation into account. It is well known that the corresponding QCD corrections can be 



^In fact, this has been the default for all modern PDF sets; see Ref. [20] though. 
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Figure 2: Representative diagrams for each of the two leading order channels. 

large. In this paper we present them for distributions of the Higgs boson. Since it has been 
shown for the inclusive cross section that the dominant SUSY effects can be approximated 
to high accuracy by an effective bbh coupling j42j . our results are directly applicable to 
the MSSM by a trivial overall rescaling. 



In Ref. |43j . a related quantity, namely the Higgs production cross section in association 
with a single tagged b quark was studied in the 5FS. Tagging a b may be useful for 
measuring the bbh Yukawa coupling yf,, for example. In this paper we consider Higgs 
production without necessarily tagging a final state b. That is, we consider the bb initial 
state, and its associated sub-channels, as a contribution to the inclusive Higgs+jet cross 
section. We present the pT and y distributions, and study the scale dependence of the cross 
section and its component channels. These results are valid to NLO in QCD perturbation 
theory. They will be presented in Section [21 

Using the knowledge of the total inclusive cross section at NNLO [28] , we can then use our 
result for NLO H+jet production in order to derive the NNLO cross section with upper 
cuts on pt- This will be described in Section [3l 

2 Next-to-leading order cross section for <t(65 — >^ /i+jet) 

The generic leading order diagrams to Higgs plus jet production are shown in Fig.[2j At 
NLO each of these receives virtual corrections, and in addition we must include the real 
emission contributions, which induce also other initial states. The full list of processes 
at NLO is: bb — > gH and bg — ?> bH at one loop, bb — > ggH, bb — >• bbH, bb — > qqH, 
gb — )■ gbH, bq — > bqH, qq — > bbH, gg — > bbH, bb — )• bbH at tree level, where q denotes 
any of the light quarks u,d,s,c. It is understood that the charge conjugated processes 
must also be included. Formally, the virtual contributions include diagrams where the 
Higgs boson is radiated off a closed bottom or top quark loop. The former lead to terms 
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Vb ' "^bZ-^H' however, which is neglected throughout our calculation, and in the spirit 



of Refs. |28p43j. we discard the latter which are proportional to the top Yukawa coupling. 
They are separately finite and gauge invariant and could simply be added to our results, 
once the ratio of top and bottom Yukawa coupling is known. 

The NLO calculation of a process as the one considered here is by now standard. We 
apply the dipole subtraction method [H] in order to cancel the infra-red poles between 
the virtual and the real radiation contributions in the bb and bg processes. Introducing the 
a-parameter for restricting the dipole phase space |45p46j not only improves the numerical 
integration, but also serves as a welcome check through the requirement of a-independence 
of the final result. Furthermore, our result for the virtual corrections agrees with the result 
of Ref. [l3]. The leading logarithmic behaviour at small px can be checked numerically 
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Figure 3: (a) Higgs transverse momentum distribution at LO (dashed) and NLO (solid); 
(b) corresponding K-factor. 



against the resummed expression of Ref. [U] . The most important check, however, is the 
numerical comparison to a fully analytic evaluation of the pT distribution to be published 
elsewhere 



In order to avoid the infra-red divergence at low Higgs transverse momenta pT, we cut 
contributions from pT < 30 GeV in this section. For our numerical analysis we use the 
following set of input parameters. The PDFs are taken from the 1VISTW2008 set [l9], and 
the QCD coupling is accordingly set to as{Mz) = 0.13939 at LO, and as{Mz) = 0.12018 at 
NLO. The bbh coupling, for which we assume the SM expression rrih/v {v = 246.22 GeV), 
is evaluated with the running bottom quark mass mb^fi^) defined in the MS scheme, 
with an input value mb{mb) = 4.2 GeV [50]. Our default value for the Higgs mass is 
Mh = 120 GeV. 

In Fig. [3] and Fig.|4] we show the transverse momentum and rapidity distributions of the 
Higgs boson, both at LO and NLO, and the corresponding X-factors 
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and E and pz are the energy and longitudinal component of the Higgs boson in the lab 
frame. The choice of the renormalization and factorization scales is given in the plots. We 
remark that the numerator/denominator in Eq. ([T]) is evaluated with NLO/LO PDFs and 
couplings. Both for the pT and the y distribution, the dependence of the K factors on 
Pt and y is very similar to what is observed for the gluon fusion channel |32fl34j : Kpj, is 
rather flat over the considered pT interval, while Ky drops mildly towards larger values of 
the rapidity. The absolute size of the corrections is significantly smaller than in the gluon 
fusion case though. Note that the K factors for the distributions cannot be immediately 
deduced from the one for the inclusive cross section for bb ^ H + X due to the strong fip 
dependence at LO, cf. Ref. [28] and Fig.[6](b). 
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Figure 4: (a) Higgs rapidity distribution at LO (dashed) and NLO (solid); (b) correspond- 
ing K-factor. Since the distribution is symmetric around y = 0, only positive values of y 
are shown. 
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Figure 5: Scale dependence of the integrated cross section for p^ > 30 GeV at Mh = 
120 GeV. (a) /iR = Mh fixed, //p varies — (b) //p = -Mh fixed, /^r varies — (c) //p = -Mh/4 
fixed, /ZR varies — (d) same as (b), but with cut on \y\ < 2. 
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Figure 6: Scale dependence of the integrated cross section for -p^ < 30 GeV at Mh 
120 GeV. (a) /iR = Mh fixed, /ip varies — (b) fip = Mh/4 fixed, /ir varies. 



In Fig. [5] we show the integrated cross section for /i+jet production with a minimum pT of 
the Higgs of pT,cut = 30 GeV: 



a{pT > PT, 



cut j 



dpT^^ — 

PT>PT,cut dpT 



(3) 



at Mh = 120 GeV as a function of (a) the factorization and (b,c) the renormalization 
scale. The factorization scale dependence is already quite small at LO and improves 
slightly at NLO. The renormalization scale dependence we show for (b) /ip = Mh and 
(c) /ip = Mh/4, respectively. For both choices, the NLO corrections improve the scale 
dependence significantly, although //p = Mh seems to lead to a more natural behaviour 
of the LO and the NLO curves. This behaviour does not change much if we restrict the 
Higgs rapidity to \y\ < 2, as shown in Fig.[5](d). 



3 NNLO cross section with Pt cut 



The knowledge of the total inclusive cross section cjtot at NNLO [28] allows us to use the 
results of the current paper to evaluate the inclusive cross section when applying a finite 
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Figure 7: Inclusive cross section with an upper cut pr,cut on the Higgs transverse momen- 
tum, see Eq. (JH), as a function ol pT, cut- The dotted, dashed, and sohd curves show the 
LO, NLO, and the NNLO result. The arrow on the right indicates the value of the NNLO 
result without cuts, cjtot- 
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Of course, PT,cnt rnust not be too small in order not to be sensitive to the region where 
large ln{pT/Mi{) terms spoil perturbative convergence. Furthermore, if a{pT < PT,cut) is 
to be evaluated with NNLO accuracy, we have to evaluate both terms on the right side of 
Eq. (g]) with NNLO PDFs. 

Fig.[6](a) shows a{pT < PT,cut) for PT.cut = 30 GeV as a function of the factorization scale 
fiF, for fiB. = Mu, varied over the rather large interval /ip € [0.1, 10]Mh. The observations 
are very similar to the fully inclusive cross section obtained without cuts |28j which is 
shown for comparison in Fig.[6](b): Perturbation theory prefers a scale significantly below 
Mh- With this choice, the renormalization scale dependence is very weak around ^r = Mjj 
already at NLO, as shown in Fig.[6](c) and Fig.[6](d). 

Finally, Fig. [7] shows the NNLO cross section a{pT < PT,cut) as a function of pT, cut- Since 
the LO only contributes at pT = 0, it is independent on pT,cut- The NNLO corrections are 
negative with respect to NLO at small PT,cut and change sign at around PT,cut ~ 50 GeV. 

4 Conclusions and Outlook 



In this paper we have presented a first study of higher order differential distributions 
for Higgs production in bottom quark annihilation. The five-flavour scheme was used to 
calculate NLO pT and y distributions of the Higgs boson. Combination with the inclusive 



NNLO total cross section allowed us to derive the inclusive cross section with upper cuts 
on the Higgs transverse momentum at NNLO. 

Concerning the choice of the factorization scale, we find that it strongly depends on the 
observable under consideration: the value fip = Mii/4: [221[2ni[2S] seems to be favoured in 
particular when the region pT = is involved, but is less motivated otherwise. 

Although we expect that the results will be phenomenologically relevant on their own, our 
approach should be useful also for an extension to a fully differential NNLO Monte Carlo 
program for the process bb ^ H + X along the lines of Ref . [51] . 
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